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Abstract 

We study the problem of homogenization for inertial particles moving 
in a time dependent random velocity field and subject to molecular diffu- 
sion. We show that, under appropriate assumptions on the velocity field, 
the large-scale, long-time behavior of the inertial particles is governed 
by an effective diffusion equation for the position variable alone. This 
is achieved by the use of a formal multiple scales expansion in the scale 
parameter. The expansion relies on the hypoellipticity of the underlying 
diffusion. An expression for the diffusivity tensor is found and various of 
its properties are studied. The results of the formal multiscale analysis 
are justified rigorously by the use of the martingale central limit theo- 
rem. Our theoretical findings are supported by numerical investigations 
where we study the parametric dependence of the effective diffusivity on 
the various non-dimensional parameters of the problem. 



1 Introduction 

Inertial particles play an important role in various applications in science and 
engineering. Examples include planet formation, particle aggregation in rotating 
flows, atmosphere/ocean science (in particular rain initiation [101135] ). chemical 
engineering. The prominent role that inertial particles play in various scientific 
and industrial applications has triggered many theoretical investigations, see for 
example [211 [M] and the references therein. 

The starting point for many theoretical investigations concerning inertial 
particles is Stokes' law which says that the force exerted by the fluid on the 
particle is proportional to the difference between the background fluid velocity 
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and the particle velocity: 

F,it)^ (y{xit),t)-xy (1) 

Various extensions of this basic model have been considered in the literature, in 
particular by Maxey and collaborators [HI [22l US HH [Ml . In this work we 
will restrict ourselves to the analysis of particles subject to a force of the form 
([T|). together with additional molecular bombardment. 

In principle the fluid velocity v{x,t) satisfies cither the Euler or the Navier 
Stokes equations and it is obtained through direct numerical simulations (DNS). 
The solution of a Newtonian particle governed by force law ([!]) coupled to ei- 
ther the Euler or Navier-Stokes equations is analytically difficult to study and 
computationally expensive. It is hence useful to consider v{x,t) in ([T]) to be a 
given random field v{x,t) which mimics some of the features of velocity fields 
obtained from DNS; one can consider, for example, random fields whose energy 
spectrum is consistent with that of velocity fields obtained from DNS. The qual- 
itative study of Newtonian particles governed by ([1]) for given (random) velocity 
fields is very similar to the theory of turbulent diffusion [19j . which has been 
primarily developed in the case r = 0. However, relatively little is known about 
the properties of solutions in the inertial case. It is important, therefore, to 
consider simplified models for the velocity field v{x,t) which render the parti- 
cle dynamics amenable to rigorous mathematical analysis and careful numerical 
investigations. 

A model for the motion of inertial particles in two dimensions was introduced 
in [37] and analyzed in a series of papers [H [13 [H mUHl [311 [35] . This model 
consists of motion in a force field comprised of a contribution from Stokes' 
law together with molecular bombardment; the velocity field is a Gaussian, 
Markovian, divergence-free random field. This gives the equations 

Ti(t) = v{x{t), t) - x{t) + cr^(t), (2a) 

v{x,t) ^V^(j){x,t), (2b) 

^ = A0+v/QC(a^,t). (2c) 

The parameter r is the Stokes number, which is a non-dimensional measure of 
the particle inertia (essentially it is the particle relaxation time). The molecular 
diffusion coefficient is given by a > , ^(i) is white noise in R^, C(3;,t) is 
space-time Gaussian white noise and A, Q are appropriate positive, self-adjoint 
operators. Gaussian velocity fields of the form Pc)) have been considered by 
various authors in the past, in particular in the context of passive advection. 
See for example 2J and the references therein. The usefulness of random velocity 
fields of the form ([^ in simulations is that, by choosing the operators A and 
Q appropriately, we can generate random velocity fields with a given energy 
spectrum, thus creating caricatures of realistic turbulent flows. Generalizations 
to arbitrary dimension d are also possible 
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Various qualitative properties of the system ^ have been studied, such 
as existence and uniqueness of solutions and existence of a random attractor 
|36| . Furthermore, various limits of physical interest have been studied: rapid 
decorrelation in time (Kraichnan) limits [T71 [521 HZj ■ 

Diffusive scaling limits (homogenization) were studied for time independent, 
periodic in space velocity fields; thus the model is obtained from (Pa|) with 
v{x,t) = v{x) only [12l HH [31] . In these three papers it was shown that the 
rescaled process 

x'{t) := ex{t/e^) (3) 

converges in distribution in the limit as e —s- to a Brownian motion with a 
nonnegative definite effective diffusivity /C. Various properties of the effective 
diffusivity , in particular, the dependence of /C on the parameters of the prob- 
lem T, a were studied by means of formal asymptotics and extensive numerical 
simulations. The purpose of this paper is to carry out a similar analysis for the 
model problem ((2|) where the velocity field is time-dependent. 

When T — 0, i.e. the particle inertia is negligible, the equation of motion 
(|2a|) becomes 

x{t)^v{xit),t)+am- (4) 

This equation has been studied extensively in the literature [111119) . The homog- 
enization problem for ^ with velocity fields of the form ([2b | . Pc)) was studied 
in [5]. There it was shown that the rescaled process ([3]), with x{t) being the 
solution of ((4|) and v{x,t) being a finite dimensional truncation of solutions 
to (j2bp . ([2c|) . converges in distribution to a Brownian motion with a positive 
definite covariance matrix, the effective diffusivity. 

In this paper we will show that a similar result holds for the inertial particles 
problem. That is, we consider the diffusive rescaling ^ for solutions to ^ with 
v{x, t) being a Galerkin truncation of ([2c| . We show, first with the aid of formal 
multiple scale expansions and then rigorously, that the rescaled processes con- 
verges to a Brownian motion and we derive a formula for the effective diffusion 
tensor. We study various properties of the effective diffusivity as well as some 
scaling limits of physical interest. Furthermore, we analyze the dependence of 
the effective diffusivity K, on the various parameters of the problem through 
numerical simulations. In particular, we show that the effective diffusivity de- 
pends on the Stokes number r in a very complicated, highly nonlinear way; this 
leads to various interesting, physically motivated, questions. 

The generator of the Markov diffusion process corresponding to ([2|) is not a 
uniformly elliptic operator, as in the case of passive tracers, but a degenerate, 
hypoelliptic operator. This renders the proof of the homogenization theorem for 
^ quite involved, since rather sophisticated tools from the spectral theory of 
hypoelliptic operators are required-see the Appendix. 

The rest of the paper is organized as follows. In Section [2| we introduce the 
exact model that we will analyze and present some of its properties. In Section 
[3|we use the method of multiple scales to derive the homogenized equation. In 
Section [4| we study simultaneously the problems of homogenization and rapid 
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decorrelation in time. In Section 5 we present the results of numerical sim- 
ulations. Section 6 is reserved for conclusions. The rigorous homogenization 
theorem is stated and proved in the Appendix. 



2 The Model 

We will study the following model for the motion of an inertial particle in M'' 

m 

Ti{t)^v{x{t),t)~x{t)+a^(t), (5) 

where t, a > and ^(i) is a standard white noise process on M'', i.e. a mean 
zero generalized Gaussian process with 

The velocity field v{x,t) is of the form 

vix,t) ^ F{x)^l{t), (6) 

where for each fixed x, F{x) : M" ^ M'^ is an n x d matrix smooth and periodic 
as a function of x, and is a stationary generalized Ornstein-Uhlenbeck 
process on M": 

m = -]A^,{t) + iVACW- (7) 

Here is a standard Gaussian white noise process on M", which is independent 
from ^{t), S > and A, A are positive definite n x n matrices. The parameter 
S controls the correlation time of the Ornstein-Uhlenbeck process ^(i). We 
remark that one can construct a velocity field (O through a finite dimensional 
truncation of ((2c)) . Notice however that we do not assume that the velocity field 
v{x, t) is incompressible as such an assumption is not needed for the analysis. We 
will, however, restrict ourselves to incompressible velocity fields when studying 
the problem numerically in Section [5l as this case is physically interesting. 

It is sometimes more convenient for the subsequent analysis to consider the 
rescaled OU process r]{t) = ^/5n{t). Written in terms of r/(t), the equations that 
govern the motion of inertial particles become 

, F(x(t))Ti(t) 
Tx{t) = ^ ^Jl'^' - x{t) + a^{t), (8a) 

Mt) = -]Arj{t) + y|c(0, (8b) 

The velocity field that appears in ([5a|) is a mean zero stationary Gaussian 
random field with correlation time 6. It is possible to show [26' that in the limit 
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as (5 ^ (the rapid decorrelation in time limit) the solution of ([5a|) converges 
pathwise to the solution of 



TX 



x{t) = F{x{t))A-'^^/KC,{t) - x{t) + crC(t). 



(9) 



The Kraichnan-like velocity field 



v{x,t) = F{x)A 



is mean zero, Gaussian, and delta-correlated in time. We will refer to ([S]) as 
the colored velocity field model and to ([9]) as the white velocity field model. 

In this paper we will be mostly concerned with the diffusive limit of solutions 
to ([S]) and ^ . That is, we will consider the rescaled process ([3]) and study the 
limit as e ^ 0. A natural question is whether the homogenization (e 0) and 
rapid decorrelation in time [6 0) limits commute. We answer this question in 
the affirmative through using formal asymptotics as well as through numerical 
investigations. 



3 Multiple Scales Expansion for Effective Diffu- 
sivities 



In this section we will derive the homogenized equation which describes the 
motion of inertial particles at large length and time scales for both the colored 
and white noise velocity fields. The derivation of the homogenized equation is 
based on multiscale/homogenization techniques [T]. We refer to [3D] for a recent 
pedagogical introduction to such methods. 

3.1 Homogenization for the colored velocity field 

We start by rescaling the equations of motion ^ according to t i— > t/e^, x i-^ 
x/e. Using the fact that for any white noise process we have that £,{ct) — -^£,{t) 
in law we obtain: 




(10b) 



(10a) 
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We now introduce two new variables 
above equations as a first order system: 



■y/rei*^ and z'^ = x'^/e and write the 



dfit) 
dz'{t) 



1 



Y{t)dt, 



1 F{z^{tW{t) 



1 



y'{t)dt, 

r]''[t)dt ' 



dt 



;V\t) dt 



(Ua) 



■dB{t), (lib) 



1 /A 



dW{t), 



(11c) 
(lid) 



with the understanding that z^ £ T*^, a;%?;' G W^, rj^ e K" and with £_{t) = 
B{t), ((t) = W{t), Bit) and W{t) being standard Brownian motions on and 
K", respectively. The SDEs pT|) clearly exhibit the two time scales 0{t) (for x*^) 
and 0{e^) (for {y'^ , z*^ , rj'^)) . Our purpose is to homogenize over the fast variables 
(j/"^, z'^ ^rf) to obtain a closed equation which governs the evolution of a;(t), and 
is valid for e ^ 1. 



3.1.1 Multiscale expansion 

Let Xf'^'''''' {x^[t),y\t),z''{t),TfAt)} denote the solution of (HH) starting at 
{x, ?/, z, ^\ and let / : M'' X M'' X T'' X 1-^ M be a smooth function. Then 
the observable u^ix, y, z, ry, t) = E/(Xf '^'^''') satisfies the backward Kolmogorov 
equation associated with the rescaled process: 



'dt 



Here: 



Co 
Cx 

Cou,n 
Cou,y 



1 



IT 

A 



2/ ■ Vj; ) li 

11 ■ V,, 



1 



1 



1 F{z)^ 



Y6 ^ 



J/ 



a2 ^ 
2 ^ 



1 



£i u% with M^|t=o = /• 



^2/ • Vj,, 



F{z)t] 

~7r 



2 ■ VjjV,;, 



0(7,?;, 



^yWy + —Ay. 



(12) 

(13a) 
(13b) 
(13c) 
(13d) 



We use the subscript OU on the last two operators to emphasize that they are 
the generators of Ornstein-Uhlcnbeck processes in rj and y respectively. The 
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operator £o is the generator of the Markov process {z{t), y{t),ri{t)} E T'' x R'^ x 



dz{t) = -^y{t)dt, (14a) 
= J-Ei^^^^!I^dt--y{t)dt+^dB{t), (14b) 
dj]{t) = -jr]it)dt + ^dW{t). (14c) 



In the Appendix we prove that the Markov process {z{t),y(t), rjlt)} given by ((T] 
is geometrically ergodic. Hence, there exists a unique invariant measure with 
smooth density p{z,y,ri) which is the solution of the stationary Fokker-Planck 
equation 

C*oP=-^iy V,p + F{z)r, ■ Vyp) + \c*o,j^yP + \c*ou,,P = 0. (15) 



Here 



^^*ou,v--=^y{y) + '^^y and £^y„-:=V„-(A77-) + |:V„V^- 



are the formal i^-adjoints of Cou.y f^nd Cou,t], respectively. The invariant 
density p(z, y, rj) is a periodic function of z and decays rapidly as ||z|| , \\y\\ — > oo. 

In the Appendix we also prove that the operator Cq (equipped with boundary 
conditions described above) has compact resolvent in the appropriate function 
space. Consequently, Fredholm theory applies: the null space of the generator 
Cq is one dimensional and consists of constants in z^y^rj. (See the Appendix). 
Moreover the equation Cof = g, has a unique (up to constants) solution if and 
only if 

(3>„ / 9{y,z,ri)p{y,z,ri)dX = Q 



where X := T'' x R'* x M" and dX := dydzdi]. We assume that that the average 
of the velocity with respect to the invariant density p vanishes: 

{F{z)rj)^ = 0. (16) 

This is natural because it removes any effective drift making a purely diffusive 
scaling natural. The identity J^y(^CQp{y, z,ri)^ dX = implies that 

{^oy)p{y,z,il)dX = 0. 



y p{y,z,v)dx = 



Consequently: 

1 F{z)r] 1 

7^ Vs 
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Thus, the centering condition is equivalent to: 

{y)p - 0. (17) 

Let us now proceed with the derivation of the homogenized equation. We 
look for a solution of in the form of a power series in e: 

u" = uo{x, y, z, 77, t) + eui{x, y, z, 77, i) + £^^2(0;, y,z,r],t) + ... (18) 

with Uj = Ui(a;, y, z, rj, t), i — 1,2, . . . being 1-periodic in z. We substitute (|18p 
into (jl2p and obtain the following sequence of equations: 

-£0^0 = 0, (19a) 

-CqUi = Ciuo, (19b) 

r ^ duo X 

-£oti2 = £iwi , (19c) 

at 

From (|19ap . we deduce that the first term in the expansion is independent of 
the fast variables z, y, 77, i.e. uq = u{x,t). Now equation (|19bp becomes: 

We can solve this by using separation of variable^: 

Ml = ^{y,z,r]) ■ VxU, 

with 

-£o<^{y,z,v)^^y- (20) 

This is the cell problem which is posed on X. Assumption PT|) implies that 
the right hand side of the the above equation is centered with respect to the 
invariant measure of the fast process. Hence the equation is well posed (see the 
Appendix). The boundary conditions for this PDE are that $(y, z, rj) is periodic 
in z and that it belongs to L^(X, p{z, y, rj) dzdydrf), which implies sufficiently fast 
decay at infinity. 

We now proceed with l|19cp . We apply the solvability condition to obtain: 

duo ir \ 

-dt ^ ^^^'^^P 

= ^(y(g)$) :V^V:,M. (21) 

Thus, the backward Kolmogorov equation which governs the dynamics on large 
scales is 

3u 

— =/C:V,V,u, (22) 



■^In principle we should also add a term which is in the null space of Co- It is easy to show, 
however, that this term does not affect the homogenized equation and for simplicity we set it 
equal to 0. 
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where the effective diffusivity is given by 

/C=-L(y®ci>)^, (23) 

and where (g) denotes the tensor (or outer) product. Notice that only the sym- 
metric part of the effective diffusivity /C is relevant in the homogenized equation 
((22|) . However, the effective diffusivity K, itself is non-symmetric in general. We 
define 

sym(/C) = i(/C + /C^) 
3.1.2 The homogenization result 

Equation ([22|) is the Backward Kolmogorov equation corresponding to a pure 
Brownian motion. We have the following result: 

Result 3.1. For e ^ 1 and t = 0(1) the function x^(t) = ex{t/e^), where x{t) 
solves ([5]), is approximated by X{t) solving 

X = v/2sj/to(/C)/3. (24) 
with X(0) — x(0), and where (3 is a standard Brownian motion on . 



A theorem, justifying the formal approximation leading to this result, is 
proved in the Appendix, using the martingale central limit theorem. 

3.1.3 Properties of the effective diffusivity 

In this subsection we show that the effective diffusivity is non negative. This 
implies that the homogenization equation is well posed. We can show this by 
using the Dirichlet form (Theorem 6.12 in [3T)J) which shows that, for every 
sufficiently smooth f{z,y,rj), 

f{-Cof)pdX = |- (|V,/p)^ + - ((V,/)^A(V,/))^ . (25) 

Now let $ be the solution of the Poisson equation (I20|) and define := a-$, 
where a S M'^ is an arbitrary unit vector. The scalar field satisfies the Poisson 
equation 

-£o0" = ■ y. 
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We combine with (|25p . to calculate 



i.lCa — —= [ {a ■ y){a ■ ^)pdX 
Jx 



= ^(|V,</'"P)^ + ^((V„0"fA(V,r)), (26) 
> 0, 

since A is a positive definite matrix. Thus the following results holds: 

Result 3.2. The effective diffusivity matrix K, given by (j23p is positive semi- 
definite and thus the limiting backward Kolmogorov equation (|22p is well posed. 

Remark 3.1. It is not entirely straightforward to check whether the centering 
condition (fTB|) or, equivalently, ([T7)) is satisfied or not, as we don't have a formula 
for the invariant measure of the fast process-we only know that it exists. It is 
possible, however, to identify some general classes of flows v{x,t) which satisfy 
p6|) by using symmetry arguments. Consider for example the case of a parity 
invariant flow, i.e a flow satisfying the condition 

F{-z) = -F{z). (27) 

It follows from ([27]) and (fT5|) that the invariant density satisfies 

p{y,z,v) ^ p{-y,-z,v)- (28) 



It is easy to see now that ((28)) implies that (|T6|) is satisfied. Hence, the centering 
condition is satisfied for velocity fields v{x, t) that are odd functions of x. 

Remark 3.2. Even if the centering condition is not satisfied, the large-scale, 
long-time dynamics of the inertial particle is still governed by an effective Brow- 
nian motion, provided that we study the problem in a frame co-moving with 
the mean flow. Indeed, if we denote by V the mean flow, i.e. 

V^{F{z)ri),, 

then the rescaled processed a;"^ = e (x(t/e^) — Vt/e^) converges in distribution 
to a Brownian motion with covariance matrix (the effective diffusivity) given by 

lC = ^l(y-.[^v]®<^\ (29) 



with 

-Co^ = l^(y-^^v), C*oP = 0. (30) 



10 



3.2 Homogenization: the white noise velocity field 

We can use the same multiscale techniques to study the diffusive scahng of ([5aj 
for the white velocity field, equation ([5]): 



Ti{t) = F{x{t))A-^VACit) ~x + a^{t). 

After a similar calculation with the colored noise problem we find that the 
backward Kolmogorov equation which governs the dynamics on large scale is 

— = /C : V,V,u, (31) 
where the effective diffusivity is given by 

£^ (32) 

Here 

-Co^y,z)^^y, Qp^O, (33) 

with 

Co -^y^z + ^F{z)A-^AA-^F^{z) : V^V, + -Cou,y (34) 

ZT T 

The operator Cou.y is defined in ^3]) . We use the notation {)p to denote aver- 
aging over T'^ x R'' with respect to the invariant distribution p. 

3.2.1 The homogenization result 

Equation (j3ip is the backward Kolmogorov equation corresponding to a pure 
Brownian motion. Hence we have the following result: 

Result 3.3. For e << 1 and t — 0(1) the function x'^{t) = ex{t/e^), where x{t) 
solves ([ni, is approximated by X{t) solving 



X = \J2sym{iqi3, (35) 
with X{Q) = a;(0), and where (3 is a standard Brownian motion on W^. 



This result can be justified rigorously by means of the martingale central 
limit theorem, as is done for the coloured noise case in the Appendix. 
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3.2.2 Properties of the effective diffusivity 

As in the case of the colored velocity field, the covariance matrix of the limiting 
Brownian motion is nonnegative define. Indeed, let a € K'' be an arbitrary unit 
vector, define 0" := a • $ and use the Dirichlet form (Theorem 6.12 in [50] ) 

aXa = '^(\Vyr?)^+^{{yyrfF{z)A-^KA-^F{zf{Vyi-))^ 
> 0. 

Thus we have the following result: 

Result 3.4. The effective diffusivity matrix K, given by (j32p is positive semi- 
definite and thus the associated backward Kolmogorov equation is well posed. 

Remark 3.3. An important observation is that the centering condition for the 
white noise problem is always satisfied. Indeed, let X := T"* x R'*, dX := dzdy 
and use the identity 

J^yQpiy,z)dx = 

together with integrations by parts to deduce that 

yp{y^ z)dx = 0. 

This suiiices for solvability of ([55)1 . Hence, the long-time, large-scale behavior 
of solutions to ^ is always diffusive. This is to be contrasted with the case of 
the colored velocity field, where an additional condition, equation (fT6|) . has to 
be imposed to ensure diffusive large scale dynamics. 

4 White Noise Limit for the Effective Diffusivity 

Consider the rescaled equation (fTO)) . and denote its solution by x^''^{t). It is clear 
that if we first take the limit as (5 ^ and then the limit e ^ 0, then x^''^{t) 
converges to a Brownian motion with covariance matrix K. given by eqn. p2p . 
without having to impose any centering condition. A natural question arises 
as to what happens if we interchange the order with which we take the limits 
(5, e — + 0. In this section we show that the two limits commute under the 
additional assumption that the centering condition (|16p is satisfied. In particular 
we have the following result: 

Result 4.1. Let A, A be positive definite matrices that commute and assume 
that the centering condition (|16p is satisfied. Then for 5 <C 1 the effective 
diffusivity K, from the colored noise model given by (|23p admits the asymptotic 
expansion 

IC = iC + OiVS), (36) 
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where IC is given by ([5^ . 

The derivation of (|36p is based on singular perturbation analysis of the cell 
problem ([^0]) and of the stationary Fokker-Planck equation (fT5)) : see [2^f251[n] . 
We start by writing the operator Cq defined in (jl3p in the form, 



with 



>Co 7 A + -^Ai + A2 , (37) 
v<5 



A = -Coct;, M ^ ^F{z)ri -Wy, A2^ ^y■'Vz + -Cou,y■ 
VT yjT T 

Of course, £q is also of the form ([57]) with replaced with A*,j = 0, 1, 2. 

Note that ^ is the generator of a d-dimensional OU process. Hence, it has 
a one-dimensional null space which consists of constants in 77. Furthermore, the 
process generated by ^0 is geometrically ergodic and its invariant measure is 
Gaussian. Since A and A commute, the density of the unique Gaussian invariant 
measure (i.e., the solution of the equation ^/o'' — 0) is 

= iexp(-?7^A-M?7), 

where Z is the normalization constant. 

Let $ be the solution of ([^0]) . As before, we define (j)" = $ • a for an arbitrary 
unit vector a e R*^. We have that 

- Cor = ^2/ ■ a, (38) 

that 

CIp = (39) 

and that 

a-ICa^^{a-yr)p (40) 
Now we need to calculate the small 6 asymptotics of (j)" and p. 

4.1 Expansion for 0" 

We look for a solution of ([551) the form of a power series in \/S. 

0" = + + (50^ + . . . 
We substitute the above into ((38)) to obtain the following sequence of equations 

= 0, (41a) 
= Ai^, (41b) 

-Ao(b^ = Ai(j)1 + A2r + ^y ■ a. (41c) 
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From Equation (|41ap we get that (j)Q = (j)Q{z,y). In order for equation (j41bp 
to be well posed it is necessary that the right hand side of the equation is 
orthogonal to the null space of A^, i.e. that 

-^F(z)f^-Wycl,^{z,y)) =0, 

Vt /pi 

which is satisfied, since the term to be averaged is linear in 77 and is a mean 
zero Gaussian density. The solution of Equation (I41bp is 

^» = J-F{z)A-'rj ■ V,0^ + $i(z, y). 
The solvability condition for (|41c|) gives 



-rj^ (F^{z)WyWycl)^F{z)A-^)7i) + -^y ■ V,0^ + -Cou.y^^ + ^y ■ a ^ 0. 



We use the fact that 



(jf B^i) pi = -B : A-^ k VBeM''^^ 



to obtain 

- £o0o = ^2/ • a- (42) 

This is precisely the cell problem for the white noise velocity field, equation 
([55]) projected along the direction a G . Hence, the small 8 expansion of the 
solution to (l38l) is 



c^iz.y. V) - c^oi^.y) + ^ ■ V,(/.? + y)^ + 0{5), (43) 

where 4)'^{z,y) is the solution to 



4.2 Expansion for p 

We look for a solution of (f551) in the form of a power series in \/5 

p = po + VSpi +Sp2 + ... 



We substitute this expansion into ()39p and equate equal powers of S to obtain 
the following sequence of equations. 

-A*oPo = 0, (44a) 
-A*oPi = AIpo, (44b) 
-A*oP2 = AIpi+A;po. (44c) 
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From the first equation we deduce that (abusing notation) pa{z, y, rf) = Pq{z, y)p^ ■ 
The solvability condition for (j44bp is satisfied since 

-/ A\podr]= ( -^(F{z)7]-\/yPo{z,y))p''dr]^0. 
The solution of (I44bp is 

Pi{z,y,v) = P'' (^-^F{^-)A^^V ■ ^yPo{z,y) + Pi{z,y)^ . 
The solvability condition for (|44c[) is 

/ {Alpi + A;po) df^ ^ 0. (45) 

We use the expressions for pq and pi to deduce that 

/■ 11 

/ A^Poiz, y, T]) d-q = -=y ■ V ^pq{z, y) + -C*ojj pq{z, y). 

and 

/ Alpi dT] = ^F{z)A-^AA-^F^{z) : VyVyPo. 
We substitute the above expressions in (|^5)l to conclude that po{z,y) satisfies 

C*oPo = 0. (46) 

Consequently, pQ{z,y) = p{z,y), the solution of the second equation in ([33]). 
Thus, the small 6 expansion of p, the solution of ([39]) is 



p = p(z,2;)p" + %/^p" (^^F{x)A~'7j -Vypoiz^y) + p,{z,y)^ +0{5). (47) 



4.3 Proof of Result 

We have that 



a-/Ca = ^(a-#")p. (48) 



In the previous subsections we expanded and p as 

0« = 0^ + V^0°+<50^ + ••• 
p = po + + (5/92 + • • • 

and showed that 

(I^Q(.x,y,ri) = '?!>o(a:,J/), 

po{x,y,v) = p'^p{x,y), 



15 



where 0q = $-a and <&, p satisfy (|33p . Now if we substitute the series expansion 
in (l48l) we obtain: 



Since 0q does not depend on rj we can integrate over the rj variable and obtain 

a-ICa = -^{a-y(t>^)p^,,y^ + 0{VS). (49) 
= a-iCa + 0{V6) (50) 

Remark 4.1. Thus we have shown that for (5^1 the effective diffusivity of 
the colored noise problem is approximately equal to that arising from the white 
noise problem up to terms of 0{d), provided the centering condition is satisfied. 
It is straightforward to show that exactly the same result holds even when the 
centering condition is not satisfied. In this case the asymptotic analysis is done 
for equations ([M)) and ([50]) : the effective drift V vanishes in the limit (5^0. 

5 Numerical Investigations 

In this section we study the dependence of the effective diffusivity ((23|) or (|32p 
on the various parameters of the problem (Stokes number, molecular diffusivity 
etc.) by means of numerical experiments. We study equations ([8|),([9]) in two 
dimensions with the velocity field being the Taylor-Green flow, modulated in 
time by a one dimensional OU process. The equations of motion for the colored 
velocity field are 

Tx ^ ^i\/^iPtg{x) - X + a^, /t = -^^+^C, (51) 

with 

i'Tcix) = sin (xi) sin {X2). (52) 

Note that we consider the original equations (O , ([6]) , (O rather than the rescaled 
version (l8al) .The white noise model is 



a 



-1 



XW^^Ptg{x)C - i + <^^- (53) 



Here ^ are independent Gaussian white noise processes in dimensions 1 and 
2 respectively. 

Our aim is to study the dependence of the effective diffusivity on the parame- 
ters a, A and 6. The Taylor-Green Flow satisfies the parity invariance condition 
(|?f|) and consequently the centering condition is satisfied. Furthermore, 
the symmetry properties of the Taylor-Green flow imply that the two diagonal 
components of the effective diffusivity are equal, while the off-diagonal compo- 
nents vanish. For the rest of the section we will use the notation K, = /Cn = /C22 
and will refer to /C as the effective diffusivity. 
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(a) a = 0.001 (b) cr = 0.01 

Figure 1: Effective diffusivity as a function of r for ct ^ 1. 



We calculate the effective diffusivity using Monte Carlo simulations, rather 
than solving the Poisson equations ([20|) . ([33|) . The numerical solution of degen- 
erate Poisson equations of this form is an interesting problem which we leave for 
future study. We solve equations (|5ip and (|53p numerically for different realiza- 
tions of the noise and we compute the effective diffusivity using its Lagrangian 
definition 

hm U{x(t) - {xm ® {x{t) - {x{t)))) = O. 

where (•) denotes ensemble average over all driving Brownian motions. In prac- 
tice, of course, we approximate the ensemble average by a finite number of en- 
semble members. We solve the equations (|51l) . (|53p using the Euler-Marayama 
method for the a;-variables and the exact solution for the Ornstein-Uhlenbeck 
process. The Euler method for the colored noise problem has a order of strong 
convergence 1 since then noise is additive in this case 1 15) ; in the white noise 
case this reduces to order 1/2, since the noise is then multiplicative. We use 
3000 particles with fixed non random initial conditions. The initial velocity of 
the inertial particles is always taken to be 0. We integrate over 10000 time units 
with M = 10~3. 

5.1 The effect of r on the diffusivity 

First we investigate the dependence of the effective diffusivity on the Stokes 
number r for the Taylor-Green flow. We set the values oi \ = a — 6 = 1. 
Our results are presented in Figures [T] and O For comparison we also plot the 
diffusion coefficient of the free particle 

We observe that when cr <C 1 the effective diffusivity is several orders of 
magnitude greater than the molecular diffusivity, both for the colored and the 
white noise case. Furthermore, the dependence of /C on r is different when 
T <C 1 and r 1, with a crossover occuring for r = 0(1). On the other hand, 
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Iree parlicle 



(a) a = 0.1179 (b) a = 1.3895 

Figure 2: Effective diffusivity as a function of t for cr = 0(1). 



the enhancement in the diffusivity becomes much less pronounced when a is not 
very small, and essentially dissapears as a increases, see Figure [H This is to be 
expected, of course. 

5.2 The effect of a on the diffusivity 

We fix now a = X = 6 = 1 and investigate the dependence of /C on cr for various 
values of r. Our results are presented in Figures [3] and [H where for comparison 
we also plot the diffusion coefficient of the free particle . 

In Figure [3] we plot the effective diffusivity of the colored noise problem in 
the case where r = 1.0 (inertial particles) and r = (passive tracers). In both 
cases the effective diffusivity is enhanced in comparison with the one of the 
free particle problem. However, the existence of inertia enhances further the 
diffusivity. This phenomenon has been observed before |28| in the case where 
the velocity field used was again the Taylor-Green velocity field but with no 
time dependence. 

In Figure |4] we plot the effective diffusivity of the white noise problem as a 
function of a in the case where t — 1.3895 (inertial particles) and t = 0(passive 
tracers). The enhancment occurs in both cases but again the existence of in- 
ertia enhances further the diffusivity. As expected when ct >> 1 the effective 

2 

diffusivities for both inertial particles and passive tracers converge to 

5.3 The effect of a and A on the effective diffusivity 

In this subesection we investigate the dependence of /C on a and A for a = 
0.1, T = 5 = 1.0. In the limit as either a ^ oo or A — > the OU processes 
converges to 0. It is expected, therefore, that in either of these two limits the 
solution of the Stokes equation converges to the solution of 

TX = -X + (T^{t), (54) 



18 











-1=0 












<?I1 



























10"^ 10"' 10° 10' 10^ 



C7 

Figure 3: Effective diffusivity as a function of a for the colored noise problem 
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(7 

Figure 4: Effective diffusivity as a function of a for the white noise problem 
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10"= 10"' 10° 10' 10= 

a, X 

Figure 5: Effective difFusivity as a function of a, A 



and, consequently, in this limit the effective diffusivity is simply the molecu- 
lar diffusion coefficient. This result can be derived using techniques from e.g. 
(Chapter 9 in [14]). On the other hand, when either a ^ or A ^ oo, the OU 
process dominates the behavior of solutions to the Stokes equation and, con- 
sequently, the effective difFusivity is controlled by the OU process. The above 
intuition is supported by the numerical experiments presented in Figure [5l In 

2 

particular, the effective diffusivity converges to ^ when either a becomes large 
or A becomes small, and becomes unbounded in the opposite limits. 

5.4 The effect of S on the diffusivity 

In this section we study the effect of S in the effective difFusivity of the colored 
noise problem. Our results are plotted in Figure [H The values of a, A are set 
equal to 1, while t — 1.3895 and a — 0.3162. 

We expect that as <5 ^ the colored noise problem should approach the 
white noise problem. This is what we see in Figure [6l since when (5 is of 0(1) 
the value of the effective diffusivity for the colored noise problem is almost the 
same as the white noise one. The rate at which the effective diffusivity for the 
colored noise problem converge to the one for the white noise problem depends 
on the values of r, a. Indeed, as we have already seen in subsection 5.1 for small 
values of r and a there is a significant difference between the values for the two 
diffusivities when 5 ~ 0(1). 

6 Conclusions 

The problem of homogenization for inertial particles moving in a time dependent 
random velocity field was studied in this paper. It was shown, by means of 
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10"' 10"' 10° 10' 10^ 

Figure 6: Effective diffusivity as a function of | 



formal multiscale expansions as well as rigorous mathematical analysis, that 
the long-time, large-scale behavior of the particles is governed by an effective 
Brownian motion. The covariance of the limiting Brownian motion can be 
expressed in terms of the solution of an appropriate Poisson equation. 

The combined homogenization/rapid decorrelation in time for the velocity 
field limit was also studied. It was shown that the two limits commute. 

Our theoretical findings were augmented by numerical experiments in which 
the dependence of the effective diffusivity on the various parameters of the 
problem was investigated. Furthermore, various limits of physical interest-such 
as tJ — > 0, r ^ etc. -where studied. The results of our numerical experiments 
suggest that the effective diffusivity depends on the various parameters of the 
problem in a very complicated, highly nontrivial way. 

There are still many questions that remain open. We list some of them. 

• Rigorous study of the dependence of the effective diffusivity on the various 
parameters of the problem. This problem has been studied quite exten- 
sively in the context of passive tracers. Apart for this being an interesting 
problem for the point of view of the physics of the problem, it also leads to 
some very interesting issues related to the spectral theory of degenerate, 
nonsymmetric second order elliptic operators. 

• Numerical experiments for more complicated flows. It is expected that 
the amount of enhancement of the diffusivity will depend sensitively on 
the detailed properties of the incompressible, time dependent flow. 

• Proof of a homogcnization theorem for infinite dimensional OU processes, 
i.e. for the model ([2]). In this setting, questions such as the dependence of 
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the effective diffusivity on the energy spectrum and the regularity of the 
flow can be addressed. 



A Proof of the Homogenization Theorem 

Let x{t) : M+ i-^ M'' be the solution to the SDE 

Tx{t) = u{x{t),t) - x{t) + cr/3i(t), (55) 

where r, cr > 0, [3i{t) is a standard Brownian motion on . Furthermore the 
field v{x, t) -.W^ y.m+ ^ is given by 

u{x,t) = F{x)n{t). 

Here, for each fixed a;, F{x) £ M'*^" and, furthermore, F{x) is smooth and 
period 1 as a function of x. Also /i(t) : R"*" ^ M" is the solution of 



^i{t) = -6-^A^i{t) + J- VA/32(0- (56) 

where (32 (t) is a standard Brownian motion on K" and A, A are n x n positive 
definite matrices. Our goal is to prove that the rescaled process 

x'{t):^ex(^-L^ (57) 

converges weakly to a Brownian motion with variance given by (|23p . We rewrite 
([55]) , ([55]) as a system of first order SDEs 

X = -^y, (58a) 

y = J-Fix)t,--y+^(ii, (58b) 

/i = -y/i +^/^2. (58c) 

This is a Markov process for {x{t),y{t), ^(t)) on M'' x M'' x M". We let z{t) 
denote the function x{t)/Z'^ so that z(i) e T'^ = R'^/Z''. Since F is 1-pcriodic 
we may view {z(t), y{t), n{t)) as a Markov process on T"* x M'' x M". 

Theorem A.l. Let {x{t), y{t), /i(t)} he the Markov process defined through 
the solution of (|58|) . where A — I and A = XI, A > 0, cr > and assume that 
the process {z{t),y{t),ri{t)} is stationary. Assume that the vector field F(x)fi 
has zero expectation with respect to the invariant measure p{z,y, fi)dzdydfi of 
the Markov process {z{t),y(t), ii(t)}. Then the rescaled process x'^{t) converges 
weakly to a Brownian motion with covariance matrix 2sym(IC) where 

/C = / {-Co<P)(»^p{z,y,fi)dzdydfj.. (59) 

JT'^xR^'xR" 
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Here ^{z,y,fi) S L^(T'' x R'^ x , p{z,y, y)dzdyd^i\W^) is the unique-up to 
additive constants -solution of the Poisson equation 

- Co<P = -^y. (60) 

Remark A.l. The assumptions on A and A are made merely for notational 
simplicity. It is straightforward to extend the proof presented below to the case 
where A, A are not diagonal matrices, provided that they are positive definite. 

Remark A. 2. In the case where the centering condition p6)) . or equivalently 
(fT7|) . is not satisfied, then to leading order the particles move ballistically with 
an effective velocity V = {F{z)ij)p. A central limit theorem of the form of 
Theorem I A. II provides us with information on the fluctuations around the mean 
deterministic motion. See also Remark 13.21 

Remark A. 3. It is not necessary to assume that the process is started in its 
stationary distribution as it will approach this distribution exponentially fast. 
Indeed, as we prove in Proposition lA.il below, the fast process is geometrically 
ergodic. This implies that for every function i/; : T'* x M'' x R" i-^ R which does 
not grow too fast at infinity there exist constants C, 5 such that 

E {'^{z{t), y{t),fi{t))) - / -0(2, y, p)p{dzdydp) 

JT^'xR'^xR" 

where E denotes expectation with respect to the law of the process {z{t), y{t), ri{t)} 
and p{z,y, fj,)dzdydfi the unique invariant measure. We make the stationarity 
assumption to avoid some technical difficulties. 

As is usually the case with theorems of the form (jA.l|) .e.g. [H [25l [TBI ITS] , 
the proof (jA.l[) is based on the cental limit theorem for additive functionals 
of Markov processes: we apply the Ito formula to the solution of the Poisson 
equation (^0)) to decompose the rescaled process ^ into a martingale part and 
a remainder; we then employ the martingale central limit theorem [51 Ch. 7] 
to prove a central limit theorem for the martingale part and we show that the 
remainder becomes negligible in the limit as e — > 0. In order to obtain these two 
results we need to show that the fast process is ergodic and that the solution 
of the Poisson equation ([^ exists and is unique in an appropriate class of 
functions, and that it satisfies certain a priori estimates. In order to prove that 
the fast process is ergodic in a sufficiently strong sense we use results from the 
ergodic theory of hypoelliptic diffusions ^U\. In order to obtain the necessary 
estimates on the solution of the Poisson equation (|20p we use results on the 
spectral theory of hypoelliptic operators [51 [51 [71 [13] ■ Our overall approach is 
similar to the one developed in [T^ . 

For the proof of the homogenization we will need the following three technical 
results which we prove in Appendix [B] 

Proposition A.l. Let Cq be the operator defined in (|13p and assume that 
F{x) e C°°(T'';R") and a > 0. Then the process {z{t), y{t), p{t)} generated by 
Cq is geometrically ergodic. 



< Ce"**, (61) 
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Proposition A. 2. Assume that A — I, A ^ XI, A > and a > 0. Also let 
p{z,y,iJ,) be the invariant measure of the process generated by Cq. Then, for 
every a € (0, 2cr~^) and j3 G (0, 2A~^) there exists a function g{z, y, /i) e 5 (the 
Schwartz space of smooth functions with fast decay at infinity) such that 

1 1 1 1 1 1 1 1 2 

p{z,y,^i) ^ e-2\\v\r-2\\t^\\ g(z,y,p). (62) 

Proposition A.3. Let h e C°°(T'' x M'* x R") with D'^y ji e ^^(T'* x x 
jg^Mj g-e||y|| -cIImII dzdydp) for every multiindex a and every e > 0. Assume 
further that J h{z,y,ri) p{dz dy d/i) = 0, where p is the invariant measure of the 
process {z{t), y{t), p{t)}. Then there exists a solution f of the equation 

- Cof - h. (63) 

Moreover, for every a, f3 > 0, the function f satisfies 

/(z,2;,Ai)-etll''ll'+fll''irV>,2/,M), f^S. (64) 

Furthermore, for every a G (0,2cr~^), /3 G (0, 2A^^), / is unique (up to an 
additive constant) in L'^{T'^ xR'^x W'',e-°'\\y\\^-^\\^'^^^dzdyd^l). 

Proof of Theorem \A.1[ We have already shown that the centering assumption 
on the velocity field, equation (fT6|) . is equivalent to {y)p — 0. Moreover, y 
clearly satisfies the smoothness and fast decay assumptions of Proposition IA.3I 
Proposition lA. 31 applies to each component of equation (pD|) and we can conclude 
that there exists a unique smooth vector valued function $ which solves the cell 
problem and whose components satisfy estimate (j64p . 

We apply now Ito formula to ^{x{t),y{t), p{t)) with {x,y,p) solving 
and use the fact that —Cq^ = 72/ to obtain 

x%t) = ex{t/e'^) 

= £x{0) + ^ / y{s)ds 

= ex{0) - 6 mz{t/e'),yit/e'),pit/e')) - $(^(0), y(0), /.(O)))) 

/ Wy^y{s),z{s),p{s))d(3i{s) 
Jo 

V,,<i>(y(s),z(s),M(s))d/32(s). 

=: ex{0) + Rl + M^+N^-. 

Clearly lime^o e^IE|a;(0)p = 0. Furthermore, the stationarity assumption to- 
gether with Propositions IA.2I and IA.3I imply that 

mil' <Ce^^z,y,p)\\l.<Ce\ (65) 
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Consider now the martingales M^ and . According to the martingale central 
limit theorem [9l Thm. 7.1.4], in order to prove convergence of a martingale to 
a Brownian motion, it is enough to prove convergence of its quadratic variation 
in to (T^ i; is the variance of the limiting Brownian motion. This now 
follows from propositions IA.2I and IA.3[ together with the ergodic theorem for 
additive functionals of ergodic Markov processes 02]. In particular, using (•)( 
to denote the quadratic variation of a martingale, we have that 

{M')t = e2— / Vy^{x{s),y{s),ll{s))®Vy^{x{s),y{s),^l{s))ds 
^ Jo 

2 

Similarly 

^ (VM*(a^, y, ® V^$(x, y, t in L\ 

We combine the above with equation (|26|) and use the fact that 77 — ^/S^jl and 
that A, A are diagonal matrices to conclude the proof of the theorem. □ 

Remark A. 4. With a bit of extra work one can also obtain estimates on the rate 
of convergence to the limiting Brownian motion in the Wasserstein metric, as 
was done in [12] for the case of a time independent velocity field. To accomplish 
this we need to obtain appropriate pathwise estimates on the rescaled particle 
velocity y{t/t^) and the Ornstein-Uhlenbeck process ^{t/e'^). We also need to 
introduce an additional Poisson equation of the type and to apply the Ito 
formula to its solution. The Poisson equation of type (ISS)) plays the role of 
a higher order cell problem from the theory of homogenization; see ,e.g., |3] 
for the proof of an error estimate using higher order cell problems in the PDE 
setting. The argument used in [TU Thm. 2.1] is essentially a pathwise version 
of the PDE argument. We leave the details of this quantative error bound to 
the interested reader. 



B Proof of Propositions I A. II - IA.3 



In this section we prove that the operator Cq generates a geometrically er- 
godic Markov process. This means that there exists a unique invariant measure 
p{dzdyd^) of the process which has a smooth density y, fj,) with respect to 
Lebesgue measure on T'' x K*^ x R"; and that, furthermore, estimate ((6T|) holds. 
In addition, we prove some regularity properties of the invariant density and 
existence and uniqueness of solutions together with a priori estimates for the 
Poisson equation ((20|) . The proof of Proposition lA.ll follows the lines of [20] . 
The proof of Propositions IA.2I and IA.3I is based on results from [H] . 

The proof of this Proposition lA.ll is based upon three lemmas. In the first 
lemma we show that the transition probability Pt has a smooth density pt with 
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respect to Lebesgue. In the second we show that pt is everywhere positive. In 
the third we show that there exists a Lyapunov function. These three lemmas 
imply that the fast process is geometrically ergodic [23 Cor. 2.8]. 

Lemma B.l. Assume that F(z) e C°°(T'^;M"). Then the Markov process 
generated by Cq has a smooth transition probability density pt- 

Remark B.l. The density pt is the solution of the evolution Fokker-Planck 
equation 

dpt 

-m ^ 

Proof. This follows from Hormander's theorem [33l Thm. 38.15]. The Markov 
process {z[t), y{t), p{t)} generated by Cq solves the SDE 



dz 



ydt, 



dy = -F{z)pdt- -ydt + -dWi 

T T T 

dp = -Apdt + VAdW2, 



with the understanding that z e T'', y G R** and p S 
drift for this system is: 



(66a) 
(66b) 
(66c) 

The Jacobian of the 



J = 



Odxd 


Idxd 


Odxn 


^DFdxd 


I J 

— -l-dxd 


~^dxn 


Onxd 


Onxd 


— An X n 



where 



(DF) 



k,l 



In order to prove that our system is hypo-elliptic we need to span R2<i+" through 
the noise vectors and their Lie commutators with the drift. We are going to 
study two cases. The first is when a ^ 0. In this case the noise provides the 
vectors: 

e^*"* = 6ij, i — d + 1, . . .2d + n. and j — I, . . . ,2d + n. 



thus we are missing d vectors in order to span ]R^'^+". We can see that we can 

j 

1 



obtain the missing vectors e^'\ i — ^, 



d in the following way: 
, i = d+\,...,2d 



thus we obtained the vectors e'*^ i — \ . . .2d + n that span 



Remark B.2. We briefly remark on the case where (7 = 0. In this case since 
there is no noise in the equation describing the motion of x and y we need to 
span M^d+n ]-^aving the following vectors: 



Si) 



2d+\,. . .2d + n. and j = 1, 



.2d- 
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We notice that: 

fO ifj = l,...,d 

( Je(*)), = <^ i^(,_2rf)(^-2d) if j = d + 1, . . . , 2d 

[^0-2d)(i~2<i)if j ^ 2d+ 1, . . . ,j = 2d + n 

So the only way not to obtain the vectors e^*' , i = c?+ 1, . . . , 2(i is if the foUowing 
equation holds: 

Fik{z) = 0, V/ e {!,.. .71} and ke{l,...k} 

If this equation does not hold we can obtain the rest of the vectors in the exact 
same way as we did in the case where a =/= 0. 

Lemma B.2. For all Z := (z, y, n) eT^ xR'^x W\ t > and open O C T'^ x 
R''xR", the transition kernel corresponding to the Markov process {z{t), x{t), fJ,{t)} 
defined in satisfies Pt{z,0) > 0. 

Proof. The proof of this result is based on a controllability argument [351 13 HO] ■ 
We start by writing (|66p compactly in the form 



dZ = Y{Z) dt + j:dW (67) 
with YiZ) = [y, ^F{z)^i-^y,-Afi]'^, W{t) = [W,it), T^zW]^ and 




The control problem associated with ([M)) is 

For any t > 0, any a e T'^ x K'' x M", and any a+ e T'* x M'* x M" we can find 
smooth U e C^{[0, i], M'^+") such that §^ is satisfied and X(0) = a, X{t) = a+ . 
To see this set R = ((X^, M^)^). The equation for M is 

We consider this equation separately since it does not involve any of the other 
state variables. Choose M to be a C°° path such that, for the given t > 0, 

M{0) ^ a3,M{t) ^ a+ (69) 

where with as, we denote the last n components of the vectors a,a~^ . Since 
VA is positive definite is invertible and is defined by substitution and will 
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be as smooth as M and (VA) ^ hence C°°. Also C/(0) can be taken as 0. The 



equation for {X^ , #^)^) is 



(PX dX , , 



where Vi is chosen that M satisfies (|68| . Now let X be a C°° path such that, 
for that given t 

X(0) \ ^ / X(t) 



where with ai,a^ we denote the first 2d elements of the vectors a, a"*" respec- 
tively. Since a is everywhere invertible, is defined by substitution and will 
be as smooth as F{X)M and - hence C°° . Also Ui{0) can be taken as 0. 
Now note that the event 

sup \\W2{s) ~ U2{s)\\ < e 

0<s<t 

occurs with positive probability since the Wiener measure of any tube is positive. 
So the Brownian motion (32 controls the rj component of equation (j66p . Again 
note that the event 

sup \\Wi{s) ~ Ui{s)\\ < e 

0<s<t 

occurs with positive probability since the Wiener measure of any such tube is 
positive. Combining now these two results is possible to deduce the required 
open set irreducibility. □ 

Lemma B.3. Let Ai be the smallest eigenvalue of A and F — sup^gi-d ||i<"(x)||. 
Then there exists a constant f3 > such that the function V{x,y,fi) = 1 + 

Co{V[x,y,^i)) < ^V{x,y,fi)+(3 

Proof. We have that V{x,y,fi) maps the state space onto [l,oo) and that 
lim||y||n^ll^oo = +00. Moreover we have 

£o{V{x,y,„)) = Fix)f,-2ry-2\\y\\^ + ^d~^^^^Af,-^,+ '-^^^^ 

trace(A) a"^ 
< ~V{x,y,^i)+ 13, 
with 13=^ + '-^^^^ + 1. □ 
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Proof of Proposition \A.1[ It follows from the above three lemmas and [5D1 
Cor. 2.8]. □ 

Using results from [H [5] we can also derive some regularity estimates for the 
invariant density. In addition, we can show that the operator Cq, the formal 
L^-adjoint of Co, has compact resolvent and, hence, Fredholm theory applies. 

Proof of Proposition \A.°A The proof of this result is similar to the proof of [12), 
Thm. 3.1], which in turn follows the lines of [El [5]. Denote by (j)t the (random) 
flow generated by the solutions to ([66]) and by Vt the semigroup defined on finite 
measures by 

{Vt^i){A)^¥.{^io^~^){A) . (70) 

By Lemma IB. II Vt maps every measure into a measure with a smooth density 
with respect to the Lebesgue measure. It can therefore be restricted to a pos- 
itivity preserving contraction semigroup on L^(T'' x R'' x M." , dz dy dfi) . The 
generator £g of Vt is the formal L^-adjoint of Lq. 

We now define an operator K on (T'' x R"* x M" , dy djj) by closing the 
operator defined on by 

K = _e7ll2'll'+7llMll"£*e-f llyll'-lllMll" 



-Ay --A, + f3[l-^]\\f,r + a[l- — ] \\y\\^ 

d 

" 2 



+ (aW-l) (y.Wy + ^)+{\p^l) L-W 



n d 

- ayF{z,fi)----. 

Note at this point that a < 2a~^ and (3 < 2A^^ is required to make the 
coefficients of and W^W^, respectively, strictly positive. 

We can rewrite the above expression in Hormander's "sum of squares" form 

as 

2d+2n 

K= J2 + ^0 , (71) 

i=l 

with 

Xi = — 7=d,i for i — 1 . . . d, 

Xi = for i = n + 1 . . . (n + d) 

Xi = ^a(^{l- )yi-n-d ior i = {n + d + I) . . .2n + d, 

X, = ^/3(^(l-^^)M»-2„-d fori = (2n + d+l)...(2n + 2d), 

y ■ V, + -) + {\P - 1) (^/i • + -j - ay • F{z,^^) - - - 



Xo 
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Since F is C°° on the torus, it can be checked that the assumptions of [51 
Thm. 5.5] are satisfied with = 1 - - Aj^ - + \\y\\^ + Combining 
this with [51 Lem. 5.6], we see that there exists a > such that, for every 7 > 0, 
there exists a positive constant C such that 

||A"+VII <C(1|A^X/1| + 1|AVI1) , (72) 

holds for every / in the Schwartz space. Clearly, the operator A^ has compact 
resolvent. This, together with (j72l) with 7 = and [51 Prop. 5.9] imply that K 
has compact resolvent. 
Notice now that 

a II l|2 II l|2 Q II ||2 I II ||2 

Thus, e^^"^" -f HmH ig fj^Q solution of the homogeneous equation 

K*e-^M'-^M\' =0. (73) 

The compactness of the resolvent of K implies that there exists a function g 
such that 

Kg = 0. 

Estimate (j72p . together with a simple approximation argument imply that 
jjA'^'gjl < 00 for every 7 > 0, and therefore g belongs to the Schwartz space. 
Furthermore, an argument given for example in |8i Prop 3.6] shows that g must 
be positive. Since one has furthermore 

£*e-^ll^^ll'-^ll^ll'.9 = , (74) 
the invariant density p satisfies estimate (|62p. □ 

The ergodicity of the fast process, together with the above proposition enable 
us to prove the following lemma. 

Lemma B.4. Let a E (0, 2cr^^), /3 G (0, 2A^^) and let K he as in the proof of 
Proposition \A.Si Then, the kernel of K is one-dimensional. 

Proof. Let g G keri^. Then, by the same arguments as above, e~^"^"^^^"'^"^.g 
is the density of an invariant signed measure for Vt- The ergodicity of Vt 
immediately implies g oc g. □ 

Now we are ready to prove estimates on the solution of the Poisson equation 

isa. 



Proof of Proposition \A.3\ By hypoellipticity, if there exists a distribution / 
such that ([55)1 holds, then / is actually a C°° function. 

We start with the proof of existence. Fix a G (0,2ct~^), /? G (0,2A^^), 
consider the operator K* defined in (|73p . and define the function 



u{z,y,n)^h(z,y,ri)e ^H'^H 
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It is clear that if there exists / such that K* f = u, then / = 6^1'''"^+^"^" V is 
a solution to ([55)1 . Consider the operator K*K. By the considerations in the 
proof of Proposition 1X21 K*K has compact resolvent. Furthermore, the kernel 
of K*K is equal to the kernel of K, which in turn by Lemma [6.41 is equal to the 
span of g. Define Ti = {g)^ and define M to be the restriction of K*K to 7i. 
Since K*K has compact resolvent, it has a spectral gap and so M is invertible. 
Furthermore, we have that f G Ti., therefore / — KM^^u solves K* f = u and 
thus leads to a solution to ([55)1 . 

Since K* satisfies a similar bound to (ffS]) and since ||A^u|| < oo for every 
7 > 0, the bound ([M)) follows as in Proposition lA.2l The uniqueness of u in the 
class of functions under consideration follows immediately from Lemma [B.41 □ 

Remark B.3. Note that the solution / of ([63| may not be unique if we allow 
for functions that grow faster than . 

References 

[1] A. Bensoussan, J.L. Lions, and G. Papanicolaou. Asymptotic analysis of 
periodic structures. North-Holland, Amsterdam, 1978. 

[2] R.A. Carmona and F. Cerou. Transport by incompressible random velocity 
fields: simulations & mathematical conjectures. In Stochastic partial dif- 
ferential equations: six perspectives^ volume 64 of Math. Surveys Monogr., 
pages 153-181. Amer. Math. Soc, Providence, RI, 1999. 

[3] R.A. Carmona and L. Xu. Homogcnization theory for time-dependent two- 
dimensional incompressible gaussian flows. The Annals of Applied Proba- 
bility, 7(l):265-279, 1997. 

[4] D. Cioranescu and P. Donato. An Introduction to Homogcnization. Oxford 
University Press, New York, 1999. 

[5] J. -P. Eckmann and M. Hairer. Non-equilibrium statistical mechanics of 
strongly anharmonic chains of oscillators. Comm. Math. Phys., 212(1):105- 
164, 2000. 

[6] J. -P. Eckmann and M. Hairer. Spectral properties of hypoelliptic operators. 
Comm. Math. Phys., 235(2):233-253, 2003. 

[7] J-P. Eckmann, C-A. Pillet, and L. Rey-Bellet. Entropy production in 
nonlinear, thermally driven Hamiltonian systems. J. Statist. Phys., 95(1- 
2):305-331, 1999. 

[8] J. -P. Eckmann, C.-A. Pillet, and L. Rey-Bellet. Non-equilibrium statistical 
mechanics of anharmonic chains coupled to two heat baths at different 
temperatures. Comm. Math. Phys., 201(3):657~697, 1999. 



31 



[9] S.N. Ethicr and T.G. Kurtz. Markov processes. Wiley Scries in Probabil- 
ity and Mathematical Statistics: Probability and Mathematical Statistics. 
John Wiley & Sons Inc., New York, 1986. 

[10] G. Falkovich, A. Fouxon, and M.G. Stepanov. Acceleration of rain initiation 
by cloud turbulence. Nature, 419:151-154, 2002. 

[11] G. Falkovich, K. Gawcdzki, and M. Vcrgassola. Particles and fields in fluid 
turbulence. Rev. Modern Phys., 73(4):913-975, 2001. 

[12] M. Haircr and G.A. Pavliotis. Periodic homogcnization for hypoelliptic 
diffusions. J. Statist. Phys., 117(1-2):261 279, 2004. 

[13] B. Helffer and F. Nier. Hypoelliptic estimates and spectral theory for Fokker- 
Planck operators and Witten Laplacians, volume 1862 of Lecture Notes in 
Mathematics. Springer- Verlag, Berlin, 2005. 

[14] W. Horsthemkc and R. Lcfcvcr. Noise-induced transitions, volume 15 of 
Springer Series in Synergetics. Springer- Verlag, Berlin, 1984. Theory and 
applications in physics, chemistry, and biology. 

[15] P.E. Kloeden and E. Platen. Numerical solution of stochastic differential 
equations, volume 23 of Applications of Mathematics (New York). Springer- 
Verlag, Berlin, 1992. 

[16] Tomasz Komorowski and Stcfano OUa. On homogcnization of time- 
dependent random flows. Probab. Theory Related Fields, 121(1):98-116, 
2001. 

[17] R. Kupferman, G. A. Pavliotis, and A. M. Stuart. Ito versus Stratonovich 
white-noise limits for systems with inertia and colored multiplicative noise. 
Phys. Rev. E (3), 70(3):036120, 9, 2004. 

[18] C. Landim, S. 011a, and H. T. Yau. Convection-diffusion equation 
with space-time ergodic random flow. Probab. Theory Related Fields, 
112(2):203-220, 1998. 

[19] A.J. Majda and P.R. Kramer. Simplified models for turbulent diffusion: 
Theory, numerical modelling and physical phenomena. Physics Reports, 
314:237-574, 1999. 

[20] J.C. Mattingly and A. M. Stuart. Geometric ergodicity of some hypo- 
elliptic diffusions for particle motions. Markov Processes and Related Fields, 
8(2):199-214, 2002. 

[21] M.R. Maxey. Gravitational settling of aerosol particles in homogeneous 
turbulence and random flow fields. J. Fluid Mech., 174:441-465, 1987. 

[22] M.R. Maxcy. The motion of small spherical particles in a cellular flow fleld. 
Phys. Fluids, 30(7):1915-1928, 1987. 



32 



[23] M.R. Maxcy. On the advcction of spherical and nonspherical particles in a 
nonuniform flow. Philos. Trans. Roy. Soc. London Ser. A, 333(1631):289- 
307, 1990. 

[24] M.R. Maxey and J.J. Riley. Equation of motion for a small rigid sphere in 
a nonuniform flow. Phys. Fluids, 26:883-889, 1983. 

[25] G. C. Papanicolaou. D.W. Stroock, and S. R. S. Varadhan. Martingale 
approach to some limit theorems. In Papers from the Duke Turbulence 
Conference (Duke Univ., Durham, N.C., 1976), Paper No. 6, pages ii+120 
pp. Duke Univ. Math. Ser., Vol. III. Duke Univ., Durham, N.C., 1977. 

[26] G. A. Pavliotis and A. M. Stuart. White noise limits for inertial particles in 
a random field. Multiscale Model. SimuL, l(4):527-533 (electronic), 2003. 

[27] G. A. Pavliotis and A. M. Stuart. Analysis of white noise limits for stochas- 
tic systems with two fast relaxation times. Multiscale Model. SimuL, 4(1):1- 
35 (electronic), 2005. 

[28] G. A. Pavliotis and A. M. Stuart. Periodic homogenization for inertial 
particles. Phys. D, 204(3-4): 161-187, 2005. 

[29] G.A. Pavliotis. A multiscale approach to Brownian motors. Phys. Lett. A, 
344:331-345, 2005. 

[30] G.A. Pavliotis and A.M. Stuart. An introduction to Multiscale Methods. 

Lecture Notes, 2006. 

[31] G.A. Pavliotis, A.M. Stuart, and L. Band. Monte Carlo studies of effective 
diffusivities for inertial particles. In Monte Carlo and quasi-Monte Carlo 
methods 2004, pages 431-441. Springer, Berlin, 2006. 

[32] D. Revuz and M. Yor. Continuous martingales and Brownian motion, vol- 
ume 293 of Grundlehren der Mathematischen Wissenschaften [Fundamen- 
tal Principles of Mathematical Sciences]. Springer- Verlag, Berlin, third 
edition, 1999. 

[33] L. C. G. Rogers and David Williams. Diffusions, Markov processes, and 
martingales. Vol. 2. Cambridge Mathematical Library. Cambridge Univer- 
sity Press, Cambridge, 2000. 

[34] J. Rubin, C. K. R. T. Jones, and M. Maxey. Settling and asymptotic motion 
of aerosol particles in a cellular flow field. J. Nonlinear Sci., 5(4):337-358, 
1995. 

[35] R. A. Shaw. Particle-turbulence interactions in atmosphere clouds. Annu. 
Rev. Fluid Meek., 35:183-227, 2003. 

[36] H. Sigurgeirsson and A. M. Stuart. Inertial particles in a random field. 
Stoch. Dyn., 2(2):295-310, 2002. 



33 



[37] H. Sigurgeirsson and A. M. Stuart. A model for preferential concentration. 
Phys. Fluids, 14(12):4352 4361, 2002. 

[38] D.W. Stroock and S.R.S. Varadhan. On the support of diffusion pro- 
cesses with applications to the strong maximum principle. In Proceedings 
of the Sixth Berkeley Symposium on Mathematical Statistics and Proba- 
bility (Univ. California, Berkeley, Calif, 1970/1971), Vol. Ill: Probability 
theory, pages 333-359, Berkeley, Calif., 1972. Univ. California Press. 

[39] L.P. Wang, M.R. Maxey, T.D. Burton, and D.E. Stock. Chaotic dynamics 
of particle dispersion in fluids. Phys. Fluids A, 4(8): 1789-1804, 1992. 



34 



